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The Casimir and electromagnetic interactions between objects at short separations are strongly 
influenced by the local geometry near the point of closest approach. In this paper we demon¬ 
strate that the assumptions underlying common statistical analysis of roughness may not hold in 
experiments using micro-spheres as interacting objects. Based on an extensive experimental and 
numerical analysis of the surface topology of the widely used colloidal particle types 4310A and 
4320A, we show that the actual variation in the local surface curvature may give rise to large 
uncertainties in the comparison of experimental data to theories. 
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I Introduction 

Over the past two decades hydrodynamic |jQ and 
Casimir EHU (or van der Waals) interactions at short 
surface separations have received increased attention. 
Numerous measurements using a wide variety of exper¬ 
imental setups have been conducted in order to better 
understand these effects. In particular, the technical 
maturity and availability of atomic force microscopes 
(AFM) has made these instruments popular for the 
measurement of forces at short surface separations. As 
both the hydrodynamic and the Casimir forces scale 
with distance as well as with the area of the opposing 
surfaces, one would ideally use a parallel plate config¬ 
uration in order to maximize the force. Due to the 
technical difficulty of maintaining parallelism, how¬ 
ever, measurements are mostly performed in geome¬ 
tries with higher degrees of symmetry such as a sphere 
opposing a flat plate or orthogonal cylinders. 

When comparing experimental results to theoretical 
predictions, the imperfections of the experimental sur¬ 
faces have to be taken into account. For hydrody¬ 
namics, a wide spectrum of influences from tribological 
parameters, such as wettability [5], surface contamina¬ 
tion mm, trapping of nano bubbles mm, and rough¬ 
ness [IMS] have been investigated (for a recent review 
see CD- In the field of Casimir physics, dielectric prop¬ 
erties men, the thickness of surface layers [19], lo¬ 
cal variations in the Fermi potential (patches) [2014221 , 
and roughness [231426] have been shown to influence 
the measured forces significantly. Stochastic irregular¬ 
ities also seem to play an important role for capillary 
forces [27] and the mechanical characteristics of micro¬ 
electromechanical switches [281 . 

Being short ranged, hydrodynamic slip effects as well 
as the Casimir force crucially depend on the geomet¬ 
rical properties of the interacting objects. In most 
studies, geometry is accounted for by considering the 
global radius of curvature and the roughness ampli¬ 
tude [29]. While these two parameters represent the 
extreme scales of surface geometries, the intermediate 
scale, namely local variations in the curvature (some¬ 


times referred to as ‘waviness’), are mostly neglected 
in studies utilizing spherical (or cylindrical) probes. 
However, theoretical studies [80] have shown that pre¬ 
cision Casimir experiments performed with lenses of 
centimeter-size radius are highly sensitive for non¬ 
sphericity. Other authors report variations of several 
tens of percents in the hydrodynamic drainage force 
between a sphere and a plate for local dents m, large- 
scale asperities [32] . and generic stochastic variations 
of the surface topology [581 . 

It is the aim of the present study to demonstrate that 
the basic assumptions of commonly applied statistical 
methods for the description of roughness (and general 
corrugations) of spheres may not always be met. Based 
on a series of atomic force microscope (AFM) measure¬ 
ments on commonly used colloidal probes, we derive 
global and local values for the surface curvature. At 
hand of experimental topology data we show that the 
Casimir force may vary by several percent in depen¬ 
dence on the position of the point of closest approach 
on the sphere - an effect which cannot be covered sta¬ 
tistically. 

The paper is organized as follows. Section [II] gives a 
theoretical and numerical assessment of the effective 
area on colloid surfaces that give the major contribu¬ 
tion to the Casimir and electrostatic interactions in 
an experiment. In Section III we describe the experi¬ 
mental methods used in the present work. A detailed 
presentation and analysis of the data on radii and vari¬ 
ations of the Casimir force in Section [TV] is followed by 
a brief summary and conclusion of our work in Sec¬ 
tion El 


II Locality of surface interactions 

Under the sole premise that the smallest curvature R 
of the interacting surfaces is much larger than the sep¬ 
aration d 0 at the point of closest approach (d 0 R), 

the Casimir force can be estimated by the Derjaguin 
approximation (also known as proximity force approx¬ 
imation (PFA)) [34] . In order to quantify the ‘locality’ 
of the interaction it is instructive to compute the size 
of the cap on an ideal sphere that contributes most 
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Figure 1 : Geometry of a corrugated sphere opposing an 
ideally flat plate (not to scale). 


to the force between the same sphere and a flat plate. 
For this purpose, we write the PFA in the form 


r c 2i r 

Fc(r c ,d) = Jdr J d iprP c ,pp(d(r,(p)) , (1) 

o o 

where d(r, ip) is the local surface separation, r c is the 
lateral radius of the cap (see Fig. [I]), ip is the azimuthal 
angle, and Pc, pp (d) stands for the Casimir pressure for 
parallel plates, as computed via Lifshitz theory [551 . 
For a perfect sphere of radius R , the distance function 
is d(r, ip) = d(r) = do + R — \/R? — r 2 , and we can 
compute numerically the ratio Fc(r c , do)/Fc(R, do). 
A parametric analysis evaluating Eqn. 0 for room 
temperature, and gold as material for all surfaces 
(where the cliclelectric function has been computed 
on the basis of tabulary data and the Drude model 
as described for example in [36]), R = 50 /jm, and 
do = 100 nm, a fraction of 95.5% of the total force 
is generated within a radius of r c = r 3a c = 4.65 pm 
around the point of closest approach. For a sphere 
with R = 100 /im, a similar estimation results in 
V 3 cr,c = 6.47 pm. Fig. i shows the dependence of 
f 3 <T,C on distance. Although being based on an ap¬ 
proximation, these results for r 3cr c clearly demon¬ 
strate that, in typical AFM Casimir experiments using 
nretal-coated colloid particles, the major contribution 
to the measurement comes from a very small fraction 
of the total surface area. We will discuss the implica¬ 
tions of these findings with respect to roughness cor¬ 
rections further in Section HV. 2 1 

The surface separation is determined in many experi¬ 
ments from a calibration based on electrostatic forces. 
Hence, it is interesting to analyze the ‘locality’ of 
the latter interaction in a similar way as described 
above for the Casimir force. For this purpose, one 
can use a modified version of the widely used theory 
by Smythe m for the electrostatic attraction between 
a plate and a sphere. According to the method of mir¬ 
ror charges, the potential </>(x) in three-dimensional 
Cartesian coordinates x = (x, y, z) can be expressed 


as, 


oo 

<M X ) = ~ A — V] 

47T£ 0 L 

j=i 


dj 


Pi 


S'? x r. ; - X 


, ( 2 ) 


where s j = ( 0,0, do + R 


1 sinh(j — 1 )a 

jsinhja J J 


dj = ~Pj = 47T£ 0 rfo (V3+hi sinh jaj 
r j = —s j , and a = cosh -1 (l + . 


The force generated by a local cap of radius r c , which 
is centered around the point of closest approach (see 
Fig. 0 is obtained from Coulomb’s law by integrat¬ 
ing the square of the electric field E = — V0(x,e? o ) 
(standing parallel to the surface normal vector n of 
the sphere, with ||n|| = 1), over the cap area, 


FE{r c ) = 27 t 


B(r 0 ) 


d0 (n 


■ E) 2 f? 2 sin 6 


(3) 


where we have used spherical coordinates (R,6,ip) 
with the relations x = R cos ip sin 0, y = R simp sin 0 , 
z = do + R{ 1 — cos0), the unit vector in z-direction 
e z with || e_ || = 1, and the definition 0(r c ) =n — 
sin -1 (r c /i?). Eqn. 0 is amenable for numerical 
evaluation and allows to determine the effective ra¬ 
dius t 3(J e for the electrostatic force from the ratio 
Fe(v c , do)/FE(R, do0 The results of this analysis are 



do [pm] 

Figure 2 : (color online) Dependence of the effective cap 
radii r 3 and r 3lT} E for Casimir and electrostatic interac¬ 
tions on distance d, for particles of radii R = 50 pm and 
R = 100 pm. 


collected in Fig. [2j At very large do, the effective radii 
contributing to both Casimir and electrostatic interac¬ 
tions converge to R. For do < 1 pm, r 3(J ^E is approx¬ 
imately a factor of 3 larger than r 3(Ti c We note that 

1. Note that the error from neglecting the upper half of the 
sphere is < 5 % for do < 3 p,m, but is irrelevant at all distances 
due to the normalization by Fe(R, do). The same computation 
with Fe ,compl (.Ri ^o) (obtained by setting 0(r c ) -A 0 in Eqn. |3j)) 
instead of FE(R,do) yields slightly smaller values of e &t 
large distance. 
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the results obtained here for ro a ,c state a non-trivial 
extension of the estimates of the effective length scale 
given in Ref. l38l . which are limited to a power-law 
description of the forces. 

Ill Experimental methods 

We have investigated the commercially available col¬ 
loidal particles 4310A and 4320A from Duke Scientific, 
which have found application in numerous experiments 
using AFM. These particles have NIST-certified diam¬ 
eters of 201±3.2 /irn^] (standard deviation 7.8 /im) and 
100 ± 1.5 (1.6) pm, respectively. 

For both types of spheres large dense arrays of parti¬ 
cles were fixed with epoxy glue on a glass slide. Small 
conglomerates of 4 x 4 colloids with optically clean 
surfaces were then selected for further investigation. 
The surface topology was measured on several spots 
per sphere by using a Veeco Nanoscope III AFM in 
closed loop tapping mode. Then, the arrays were RF 
sputtered with 5nm Cr followed by 100 nrn Au, fol¬ 
lowing the same procedures and parameters used for 
the preparation of probes for recent Casimir and hy¬ 
drodynamic experiments mmm- All spheres were 
imaged in a FEI Phenorn SEM, followed by a second 
series of measurements in the AFM. 

For the purpose of this study, a proper calibration of 
the instruments is of vital importance. Prior to the 
first AFM measurements we performed the full cali¬ 
bration procedure following the steps provided by the 
manufacturer using a commercially available NT-MDT 
TGG1 grating, which consists of a flat surface with 
well defined periodical triangular ridges. As a test 
of the calibration, and in order to give a proof for 
the feasibility to accurately measure the surface curva¬ 
ture radius, we have performed tapping mode imaging 
on stripped Corning SMF-28 fibers featuring a very 
precisely defined nominal diameter of 125 ± 0.7/tin 
with a maximum non-circularity of 0.5%. By fit¬ 
ting the resulting topology data of the cylindrical 
surfaces we obtained a value for the radius equal to 
62.61 ± 0.01 (0.39) /im. This result demonstrates the 
fitness of the applied method to accurately determine 
the curvature radius of a surface. For clarity, we have 
to note that radius values obtained from surface height 
data throughout this work are based on unweighted 


least square fits (see Section IV. 1). We perform statis¬ 
tical averaging of radii measured on the same sphere 
(or fiber), using the inverse square of the 68.3% con¬ 
fidence intervals of the parameter R of the fits as 
weights, and report the weighted mean, its uncertainty, 
and the weighted standard deviation (in brackets). 
Tapping mode scanj^] are performed along a rectan¬ 
gular raster, resulting in a fast (along the single lines 
of the raster, Xf) and a slow (orthogonal to the lines, 
x s lx/) scan direction. In the evaluation of the topol¬ 
ogy of the SMF-28 samples only data have been taken 
into account which were acquired with the fast scan¬ 
ning direction oriented orthogonally to the fiber. Scan¬ 
ning parallel to the cylinder axis (with x s being orthog¬ 


onal to the fiber) resulted in deviations from the stated 
results due to the presence of drift. We found that 
this drift was rather constant over all our measure¬ 
ments and samples, yielding a distance scaling factor 
fc = 0.70 ±0.06 between calibrations performed along 
x s and xj. respectively. A rotation of the scan direc¬ 
tion by 90 degrees without moving the sample yielded 
the same value for /. Finally, using the TGG1 grat¬ 
ing, we could completely exclude the appearance of 
drift along Xf. Therefore, we speculate that this sys¬ 
tematic effect is caused by an error in the scanner of 
our AFM. However, due to the strong reproducibility, 
consistency, and linearity, the drift can be removed 
from the data by the application of a scaling factor 


along x s (see Section IV). 


IV Analysis and Results 

Several series of analysis have been performed on all 
data. First, in Section [IV. 1| the accuracy and spread 
in the determination of sphere radii from AFM scans 
are investigated. The results are then compared to 
the outcomes of a more commonly used method using 
a SEM. In a second analysis in Section IV. 2[ we focus 
on the influence of locality on the Casimir force, and 
the local variation in the surface curvature. 


IV. 1 Radius determination via AFM 

In order to characterize the investigated sample of par¬ 
ticles statistically, and to further verify our method by 
comparison of the results to the literature [22] > the cur¬ 
vature radius of 32 spheres has been measured. For 
this purpose, AFM topology data hAFivi{x,y ) with 
lateral dimensions 15 x 15 /im 2 have been analyzed 
without the application of filters or any automatic cor¬ 
rections in order to avoid the introduction of system¬ 
atic errors. Using a standard least squares method, 
h'AFM (u y) was fit by a function 

hf{x,y) = h 0 + R- \/R 2 -(x-xo) 2 -(fDy-yo) 2 , 

(4) 


where x and y are the fixed lateral coordinates of 
the scan, R. Xq , yo , and ho are free parameters 
representing the global radius of the sphere, and 
offsets in all three Cartesian coordinates, respectively. 
fo is a scaling parameter, which accounts for the 
constant drift in the slow scanning direction (see 
discussion in Section III). Note that, because fn 
also is a free fit parameter, the possibly erroneous 
curvature along the slow scanning direction y = x s 
does not influence the final result for R. By averaging 


2. Note that NIST-certified values are based on only two mea¬ 
surements and may, thus, not be statistically representative. 

3. We would like to note that the investigation of spherical 
particles by means of standard tapping mode measurements has 
proven to be prone to approach errors, which eventually lead to 
damage on the probe surface. For future experiments, it is thus 
recommended to use the reverse imaging technique described 
in the literature m to characterize the surface of colloidal tip 
AFM probes. 
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Table 1 : Results for the sphere radius R, and the roughness a (weighted averaged mean over results from scans on 
the same sphere type, format: RMS{PV}), as obtained by SEM and AFM measurements. 


type SEM AFM before dep. AFM after dep. 

R [/zm] R [/tm] a [run] R [/zm] cr [nm] 

4310A 49.8±0.2 (0.7) 48.93 ± 0.46 (1.3) L7 {35} 47.59 ± 0.50 (1.4) 5.7 {137} 

4320A 102.3A0.5 (3.6) 100.9 ±2.8 (6.9) 10.2 {131} 99.3 ±3.0 (7.3) 13.2 {167} 


over the outcomes of all performed data fits, we obtain 
a value fn = 0.75 ± 0.05, which compares very well to 
the drift factor /c, determined in the calibration, as 
mentioned in Section IlIIl 

The root mean square (RMS) and peak-valley (PV) 
roughness have been measured from the fit residuals 
hAFM — hf. In order to obtain an independent set 
of measurements to which the results of our AFM 
analysis can be compared to, we have determined 
the radius of all colloids in the present study by 
means of circular fits to their circumference in SEM 
images. Results for the roughness and the global 
radius, as extracted from AFM data taken before 
and after the sputter coating, and from SEM data 
(only after the deposition) are given in Tab. [lj The 
radii extracted from SEM measurements generally 
comply with the NIST certification of the spheres. 
In direct comparison, the results for R from AFM 
fits before the sputter process show good agreement 
with the corresponding values from SEM data: 
^afm/^sem — 1 = 1.9 A 0.3 (2.1) % on 4310A and 
0.4 ± 0.4 (5.8) % on 4320A. After the deposition of 
gold, a slight change (of —5.5 ± 3.0 (2.2) %) can be 
observed in 7?afm for of spheres of type 4310A, while 
the results for 4310A remain constant within the error 
of determination (relative change 0.6 ± 0.5(10.5)%). 
The average roughness obtained from the same data 
increases after sputtering. This change is most 
pronounced for the (initially) relatively smooth 4310A 
colloids, for which the RMS value rises by a factor 
~ 3.4, and the peak value by a factor 4.1. With 
applied coating, the PV roughness of 4310A and 
4320A colloid species are comparable, while the 
RMS value still stays clearly lower on 4310A. These 
observations compare well with the findings of other 
investigations on the same sphere types [25] 02]. 

We would like to note that the method applied to 
extract the curvature from AFM topology data in the 
present work is different from the one of Refs. no on]. 
These authors used only the height profile along one 
single cross section of the data to derive the global 
radius while our fit to hf(x,y) is applied to the entire 
data set. 


IV.2 Locality, statistics, and the Casimir force 

As discussed in Section [TTJ due to the finite curvature 
the actual area contributing most to the Casimir and 
electrostatic forces between a sphere and a plate is very 
small. In consequence, it is clear that the interaction 
measured in an experiment crucially depends on the 


local topology of the surface. This issue has been in¬ 
vestigated theoretically and experimentally for hydro- 
dynamic forces in the case of a deformation (flattening) 
of the colloid [3Tj as well as for asperities [32] and pat¬ 
terned surfaces [ 32 ]. All studies have found variations 
of several tens of percent in the force. For the Casimir 
force, dedicated investigations 031 - 05 ] have shown a 
significant dependence of the interaction on the rela¬ 
tive lateral position of regular corrugations on a sphere 
with respect to an opposed flat surface with structures 
of the same periodicity. Also, it has been recognized 
early m that non-stochastic large-scale deviations 
from perfect sphericity have a significant influence on 
the Casimir force. However, for standard sphere-plate 
experiments reported in the literature, in the analy¬ 
sis of measurements the imperfections of the surfaces 
of spherical AFM probes are considered statistically 
only, either via perturbative approaches (see for exam¬ 
ple (24] 02] 08]), application of the PFA (36] 09] [50] , or 
combinations thereof [261151] . The underlying assump¬ 
tions defining the applicability of statistical methods 
are 0H] that: 

1. the area associated with the lateral correlation 
length £ (see below and Ref. [53]) of the rough¬ 
ness is much smaller than the effective interaction 
area: RdoTr £ 2 7t, 

2. the roughness profile is invariant with respect to 
the lateral position of the sphere. 

In typical analysis found in the literature, roughness 
parameters are determined from AFM scans, which are 
flattened in order to remove the curvature (in the case 
of a sphere) and to compensate for piezo drift (for open 
loop scans). Depending on the applied method of flat¬ 
tening and the size A of the scan area, low-frequency 
spectral components are eliminated from the profile 
h.(x,y ) (for the sake of simplicity, we drop the index 
AFM from now on). This latter process influences the 
shape of the auto-correlation function 


ACF(r) = (1/A)(^J dx 1 dy' \\h(x + x',y + y')h(x,y)\\^ , 

(5) 


from which £ can be obtained under the assumption 
of purely stochastic roughness by fitting ACF(r ) to 


the function ACF^ t {r) = T 0 exp — (r/£) 


2 H 


where 


0 < H < 1 is the ro ughness exponent, r 0 is a con¬ 
stant, r = yjx 2 + y 2 , and the notation () r indicates 
averaging over all pairs (x,y) on a circle of radius r. 
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(b) Histogram and cumulative height 


Figure 3 : (color online) Statistical data on the scan shown in Fig. [d] a) correlation functions ACF(r) as defined in 
Eqn. <§, computed from data in Fig. |4b| (flattened by subtraction of hf, dashed line), or after subtraction of a free 
polynomial of fourth order from the data in Fig. Ealsolid line), a) differential and cumulative height distribution of 


unfiltered data of the flattened AFM scan in Fig. 4b 


Fig. 3a shows an example of the function ACF(r ) ob¬ 
tained from unfiltered AFM data shown in Fig. |4] If 
the best fit of h(x, y) to a sphere according to Eqn. © 
is removed from the data (resulting in the topology 
shown in Fig. [4b|, a subsequent fit of ACF(r) results 
in £ = 2.5 /im (dashed black line in Fig. [3a] ) . De¬ 
spite the fact that the height distribution (solid line 
in Fig. 3b) almost perfectly follows a Gaussian curve, 
the spectral distribution ACF{r ) in Fig. 3a cannot be 
fit satisfactorily by the model of Eqn. As can 

be seen in the figure, a significant contribution comes 
from large-scale corrugations with wavelengths A corr 
of several /mi[^J X corr £. These contributions of low 
spatial frequency break the conditions (1.) and (2.). 
Consider for example an experiment with R = 50 /.an, 
and £ = 2.5/an. At a separation do = 100 nm roughly 
three structures of area £ 2 7t fit into a cap of size 
and Rdo/t? = 0.8. Therefore, a statistical evalua¬ 
tion of roughness might in this case not be rectified 
at do 1$ 1 M m - It is interesting to note that, if one 
removes independent polynomials of fourth order in 
both lateral directions instead of the spherical surface 
hf , the fit to ACF leads to £ = 420 nm (solid light blue 
line in Fig. 3a). The latter value compares well with 
those reported in the literature (200 nm in Ref. [50] 
and < 600nm in Ref. [56]). 


It is further instructive to perform explicit calcu¬ 
lations of the Casimir interaction between an ideally 
flat surface and the height profile h(x, y) of a sphere 
obtained from actual AFM data. We select a repre¬ 
sentative scan, which was recorded on a 4310A sphere 
after sputtering. Fig. 4a shows the topology h(x, y) [^] 


where, for the purpose of a more clear illustration, a 
running mean filter with a lateral width of 100 nm has 
been applied. For the unfiltered data, we calculate a 
statistical (stochastic) roughness correction Xr(do) to 
the Casimir force based on the established model of 


Ref. [12, which applies the PFA. For the ease of com¬ 
putation we use a BSpline-interpolation of the differen¬ 
tial height quantity w(h) of the histogram in Fig. 3b 
(based on the flattened topography in Fig. 4b), and 
write 


/*max(/i) 

Xr(d 0 ) = / dhw(h) 

J min(h) 

/*max(/i) 

with 1 = / d hw(h). 

J min(h) 


F c , PP (d 0 H) 


Fc ,pp (do) 


( 6 ) 


Results of Eqn. 0 will be taken as a reference for the 
subsequent analysis. 

From the unfiltered data of the scan shown in Fig. 4a 
we extract 36 equally distributed caps of radius r c = 
4/im and center positions ( x c ,y c ) |^J For each cap 
c, the height profile h c has been transformed such 
that the vertical 2 -axis is parallel to the normal vec¬ 
tor associated with the global spherical fit at (x c ,y c ) 
(see Fig. |5|. The Casimir force, which would act be¬ 
tween the curved caps with a corrugated surface pro¬ 
file h c (x, y) and a perfectly flat plate (in this example, 
both made of gold) at a distance do can be computed 


4. We note that the large-scale deformations, which have 
a typical height of several nm (4.2nm(RMS), at wavelengths 
> £ = 2.5 ^m), cannot be explained by piezo drift since all mea¬ 
surements have been performed in closed loop operation with 
accuracy better than 0.5 nm. Moreover we could not find any 
regular patterns or correlations between the irregularities shown 
by different scans. 

5. Note that the scan in Fig. [4] was taken from an area of 
15 x 15p,m 2 . The deviation of the aspect ratio from unity is a 
direct consequence of the scaling factor fr> (here in the direction 
of x). 

6. We would like to note that the same investigations have 
been repeated for several different spheres. In all cases, we ob¬ 
tained equal qualitative (but not the same quantitative) results 
for all points of the presented discussion. 
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(a) original topology 




(b) flattened topology 


Figure 4 : (color online) Topology of an exemplary 4310A 
sphere after sputtering, a) Actual AFM data (h afm), b) 
flattened by subtraction of the spherical fit hf. 



AFM scan x,y 


Figure 5 : Adaption of the local coordinate system of the 
cap to the tangential plane to the sphere at (a : c ,y c ). 


important aspect for the application of Eqn. ([7]) is the 
definition of the reference do for the surface separation. 
In order to demonstrate the impact of the latter pa¬ 
rameter, we elaborate the data in three different ways: 



51 51.5 52 52.5 53 53.5 54 54.5 55 

R c [pm] 


Figure 6 : (color online) Definitions and results for global 
radius reference, a) Definition of the distance reference do 
from the global spherical fit (dashed line, R), b) results for 
the normalized Casimir force Fc, c /Fc,id.{R,r c ) as a func¬ 
tion of do, computed for 36 caps of r c = 4 pm on the AFM 
scan shown in Fig. |4a| c) inverse dependence of the force 
on the local curvature radius R c resulting from higher tips 
of more convex areas. Each data point corresponds to a 
different cap. 

Electrostatic distance determination. In experi¬ 
ments where the distance between the sphere and the 
plate is determined via an electrostatic calibration, the 
surface separation do corresponds Q physically to the 


by numerical application of the PFA, 

Fc,c{do) = ~j( 7 ~ ^ ' Pc,pp (d(x n , yn)) > (7) 

c n 

where n runs over all N c points of the cap, d(x, y) = 
do — h c (x , y), and A c = is the area of the cap. An 


7. The Smythe theory for the electrostatic interaction between 
a sphere and a plate, which is used in this analysis as well as in 
the literature, does not take into account local variations in the 
solution of the Poisson equation due to inhomogeneous conduc¬ 
tivity (patch effects) and corrugations of the surface. We believe 
that this is an important point of discussion which, however, 
goes beyond the scope of this paper, and we restrain ourselves 
to the investigation of the influence of different hypothesis found 
in the literature regarding the surface separation. 
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gap between the surface of an ideal sphere, fitted to a 
cap of the effective radius r 3 <j t E, and the flat plate. In 
the analysis of the locality of electrostatic and Casimir 
forces in Section 0(Fig.|) it could be seen that t , 3 (t,e 
is on the order of 10 pm in the distance range acces¬ 
sible by an AFM experiment (50-400 nm). This in¬ 
dicates that the cap size over which the electrostatic 
interaction is integrated (averaged), is comparable to 
the one of the entire scan shown in Fig. |4j Therefore, 
we define here do between the surface of the spherical 
fit with curvature R to the entire scan (global fit), and 
the ideally flat opposing plate, as depicted in Fig. [6}r. 
However, the Casimir force is still (effectively) being 
generated within a cap of radius r 3a ^c < f 3 <j,E on the 
curved surface. When Fc,c(do) is computed for each 
spherical cap, according to Eqn. 0 one obtains the 
force curves shown in Fig. 0- Since we are interested 
in the deviation of the forces Fc, c {do) from a theoret¬ 
ical prediction Fc,id(R, i" c , do) for a cap of the same r c 
on an ideally shaped sphere of radius R , we normalize 
all given results by this latter quantity. In Fig. 0 , 
the geometric mean and standard deviation over the 
36 caps are indicated as well as the stochastic correc¬ 
tion Xr{do)- Note that, due to the lateral separation of 
1 /jm between the center positions of adjacent caps and 
the resulting overlap, the force curves are not strictly 
statistically independent. Nonetheless, the statistics 
are representative for the results in an actual exper¬ 
iment, in which the point of closest approach would 
be shifted in the same way by turning the sphere. The 
large spread of ~ ±35 % in Fc, c (do) at the shortest sur¬ 
face separation can intuitively be understood from the 
fact that the top surfaces of some caps lie closer to the 
opposing surface than others, resulting in a deviation 
between the effective surface separations for Casimir 
and electrostatic interactions. The caps indicated by 
‘cl’ and ‘c2’ (in Fig. [4]) correspond to the minimum 
and maximum force lines, respectively, in Fig. and 
represent the caps with the lowest local minimum and 
highest maximum values of h c , respectively 0 An in¬ 
teresting fact to note is that ‘cl’ is a convex area in 
Fig. [4j which indicates that I? c (cl) < R , while ‘c2’ is 
concave R c { c2) > R. The linear dependence on the 
curvature radius R, as predicted by the PFA [S3], is 
thus overruled by the influence of the distance shift, 
as shown in Fig. [6];. 


Distance determination from contact. If the 

distance is adapted to the local curvature radius R c 
(derived from a fit to the local spherical surface within 
r c around patches on the scan shown in Fig. 4a), the 
separation between the sphere and the plate has to be 
redefined as depicted in Fig. 0 - According to Ref. (38: 
this distance reference corresponds to the one that 
would be obtained experimentally by moving the sur¬ 
faces to contact, and subsequent consideration of the 
statistics of the roughness distribution measured on a 
flattened surface. We note that the local adaption of 
do to the topology of caps of the size of r 3rJt c (instead 



50 100 150 200 250 

d 0 [nm] 



R c [pm] 


Figure 7 : (color online) Definitions and results for local 
radius reference, a) Definition of the distance reference 
d'o = do — ho from the local spherical fit (dashed line, R c ), 
b) normalized Casimir force Fc, c /Fc,id{R, r c ) as a function 
of d'o, computed for 36 caps of r c = 4pm on the AFM scan 
shown in Fig. [d] c) almost linear dependence of the force 
on the local curvature radius R c . 


°f r 3 a,E ) naturally increases the variation in do but 
reflects more accurately the situation ‘experienced’ 
by the Casimir force. Accordingly, the variation in 
the F c , c (do) between different caps for this example 
(see Fig. 0 ) drops by a factor ~ 7 with respect to 
the one shown by the data in Fig. 0- As can be 
seen in Fig. [7J;, the linear dependence of Fq, c on 
the surface curvature R c is restored. In accordance 
with this finding, the minimum line in Fig. 0 now 
corresponds to l c2’ and the maximum line to ‘cl’, 
which is in agreement with the fact that the convex 
areas in Fig. 0 have a smaller local R c (and, thus, 
produce a weaker force) than concave areas. Note 
that the lateral distance between ‘cl’, and ‘c2’ is 
only 6.4 pm, corresponding to an angular shift of 


8 . Due to the circular shape of the caps, the upper left and 
lower right corners of the scan cannot be reached, and do not 
enter the calculations, thus. 
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7° on the sphere [^] Note further, that single force 
curves deviate not only in amplitude but also in the 
exponent of the distance. Therefore, the minimum 
and maximum lines do not necessarily correspond 
to actual force curves but show the min. and max. 
values of all curves at the same distance. 


Flattened topography, geometric average refer¬ 
ence. Finally, we investigate the influence of large- 
scale corrugations on the flattened topology shown in 
Fig. |4)r, obtained by removal of the curvature through 
subtraction of the global reference hf (corresponding 
to the global radius R) from the AFM data. As shown 
in Fig. [Sp, the reference do is defined by the geometric 
average of the local surface profile within ro a ,Ci accord¬ 
ing to the relation J2n= i[^o — d(x n ,y n )] = ( p*| If the 
forcUfl is computed between a perfectly smooth plate 
at a distance d 0 (see Fig. Ufo and caps at the same 
positions as before, but now on the flattened profile of 
Fig. |4b[ the spread is again reduced slightly. More im¬ 
portant, the change in the topology from curved to flat 
also leads to an up-bending of the curves at short dis¬ 
tance in Fig. [Bp which represents a qualitative change 
of the results with respect to the case of curved sur¬ 
faces. This effect is caused by the geometric (linear) 
method of statistical averaging of the values h(x,y). 
A height distribution following (according to the his¬ 
togram in Fig. 3b I a Gaussian statistic around do, will 
naturally, due to the non-geometric distance scaling of 
the force, lead to a distribution of Fc that is shifted 
towards higher values. Mathematically, the average 
distance do = d(x n , y n )/N results in an aver¬ 

age force F c = [F c , pp (d 0 )/N} J2n=i ( d o/d{x n ,y n )) 3 > 
F pp (do), where the distance dependence of the Casimir 
force Fc, PP between ideal parallel plates has been ap¬ 
proximated for the sake of simplicity by 1/d 3 . This 
effect of ‘shifting by weighting’ is pronounced much 
stronger for the case of the flattened surfaces in Fig. [8] 
than for curved surfaces. Again, the reason is the lo¬ 
cality of the interaction. While on caps with finite R c 
the central region contributes most, and only the very 
realization of the roughness at the point of closest ap¬ 
proach determines the actual force, on a (flattened) 
cap with R c —> oo all points contribute in the same 
way, and the statistical averaging enhances the force. 
Thus, there is a qualitative difference between the ac¬ 
tual dependence on roughness for flat and curved sur¬ 
faces. 

Numerical evaluation of Eqn. © results in the blue 
dotted lines in Figs. [6p [7p and ©>■ In all three fig¬ 
ures, the statistical prediction \r lies within the blue 
(innermost) band corresponding to the standard devi¬ 
ation of the 36 force curves at each distance. Nonethe¬ 
less, in the case of curved surfaces, Xr shows significant 
deviations from the mean of the force curves obtained 
in our explicit example. For the computations on flat 
surfaces, the local averages do(c) do not vary strongly, 
leading to a better agreement with \ r (for which the 


plate 

flattened 
(removal of R) 

Jd 0 

d(x,y) 

- 


R 



R c [pm] 


Figure 8 : (color online) Definitions and results for a 
topology flattened by removal of the global spherical fit R, 
and the local geometric mean do as reference, a) Definition 
of the distance reference do from the local geometric aver¬ 
age of the flattened cap (short dashed line), b) normalized 
Casimir force Fc lC /Fc,id{pp,r c ) as a function of do, com¬ 
puted for 36 flattened caps of r c = 4 pm on the AFM scan 
shown in Fig. [4] c) only slight dependence of the force on 
the local curvature radius R c after flattening of the profile. 


average distance do(,\r) is equal to the plane given 
by the global spherical fit indicated by R in Fig. |8fi). 
However, it has to be considered that, according to the 
analysis above, a purely statistical treatment may not 
be valid for curved surfaces. As has been indicated by 
other authors as well EH S 3 , one has to consider that 
in an actual experiment only one of these force curves 

9. The same analysis on other spheres lead to shifts of less than 
3 fim. Also, the caps giving the minimum and maximum forces 
for different height reference are not necessarily the same, as in 
the example in the main text. The general conclusions regarding 
the influence of the distance error and the local curvature hold. 

10. Note that the this situation of a flattened rough sphere 
opposing a perfectly smooth plate is equivalent to the model 
used to derive expressions for roughness corrections to Fc for 
the sphere-plate geometry |4811491151| . 

11. Note that Fq c is now normalized with respect to the force 
F Ct id(pPi r c ) between ideal flat circular plates of lateral radius 

r c - 
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will be effective. The curves corresponding to single 
realizations show significant variations in both ampli¬ 
tude and exponent of the distance. Unless the very 
point of contact on the spherical surface can be deter¬ 
mined with sub-/im accuracy and the topology at this 
very point is taken into account in a comparison with 
theory, there may always be a significant deviation of 
analytic and experimental results. 

The previous qualitative investigation was based on 
a single (representative) AFM scan and can therefore 
not be used to derive a quantitative estimate of the 
error due to local variation of the surface topology. In 
order to achieve the latter, we assume that the dis¬ 
tance is determined as outlined in Ref. [38] , which cor¬ 
responds to the case that the height reference is ob¬ 
tained from the local spherical fit. Then, the variation 
of the expected Casimir force scales directly with the 
curvature radius (see Fig. &)• We have performed fits 
using Eqn. ([ 4 ]) to caps of different r c for all available 
unfiltered AFM data. The results are shown as his¬ 
tograms in Fig. [9] For each sphere, the R c have been 
normalized with respect to the corresponding global 
fit radius R. Numerical values for the weighted mean 
R c and standard deviation a c of each distribution are 
given in tabulary form in the figures. Throughout, the 
agreement between R c and R is very good, which rules 
out systematic effects. 

Considering variations in the Casimir force, the most 
important parameter of this analysis is <JR{r 3a ^c) 
which, averaged over all spheres of type 4310A for 
r c = r 3<Jt c at do = 100 nm, gives a value 1.4% as 
compared to the standard deviation 1.8% of the 36 
caps with r c = 4 /ini used in the evaluation of the 
scan shown in Fig. [4j On 4320A colloids, (JR{r 3 a ,c) 
takes a similar value of 1.7%. The widths an of the 
(fitted Gaussian) distributions, naturally increase for 
smaller cap sizes (corresponding to smaller surface sep¬ 
arations) and approach, in the limit r c —>■ 0 the respec¬ 
tive value of the roughness evaluation using the full res¬ 
olution of the AFM scans. However, one has to bear 
in mind that these statistical measures do not apply in 
the comparison of force measurements taken between 
a single spot on a sphere, and a smooth plate. As dis¬ 
cussed above, an actual experiment will be influenced 
by the local realization of the geometry of the surface, 
and the resulting uncertainty may be much larger than 
(Jr. We would like to point out, that this potential 
error is neither covered by a typical statistical anal¬ 
ysis based on a stochastic distribution of roughness 
nor by consideration of the uncertainty in R deter¬ 
mined from SEM images. The resolution of the latter 
is clearly insufficient to quantify surface irregularities 
with a height of the order 10 nm. 

Finally, we note that the present investigation is fo¬ 
cused on the experimentally relevant situation of a 
micro-sphere opposing a flat plate, where at least one 
of the interacting surfaces has a finite curvature. In 
the case of parallel plates (as for example in mirco- 
electromechanical devices), the area of interaction is 


much larger. Therefore, the conditions (1) and (2) 
mentioned at the begin of Section IV. 2| are met, and 
the effect of surface corrugations can be modeled sat¬ 
isfactorily statistically. 


V Conclusion 

We have investigated the radius and surface roughness 
for colloid particles of type 4310A and 4320A, which 
are widely used in experiments on hydrodynamic and 
Casimir interactions in the sphere versus plate geome¬ 
try. Numerical analysis shows that on curved surfaces 
only a small fraction of the surface area effectively 
contributes to the measured Casimir and electrostatic 
forces. On the basis of actual AFM topology data, our 
study indicates that 

• surfaces of colloid particles show non-regular cor¬ 
rugations on intermediate length scales (1-10 ^m) 
of comparable or larger amplitude than short scale 
roughness. Flattening procedures, which are neces¬ 
sary to measure roughness distributions, may filter 
these irregularities. 

• in the presence of such corrugations, the transla¬ 
tional invariance and the statistical representativ- 
ity of the interacting ‘cap’ area may be violated at 
surface separations up to a few /im. Therefore, the 
fundamental assumptions for a statistical analysis 
may not hold. 

• forces, which are measured on different spots of 
the same surface vary in both amplitude and ex¬ 
ponent of the distance dependence. The scale of 
these variations depends on the applied reference 
for the distance. 

• the effective length scale for electrostatic interac¬ 
tions is roughly 3 times larger than the one for 
the Casimir force. Thus, the quantum-mechanical 
force ‘sees’ a much smaller area on the curved sur¬ 
face than does the electrostatic one. In the pres¬ 
ence of corrugations at intermediate length scales, 
this may lead to different offsets and definitions 
for the distance of the two interactions, which may 
result in a spread in measured Casimir forces of 
several tens of percent. 

• in the case that the distance is determined from 
the contact of the two surfaces, the error which 
would appear due to an electrostatic calibration 
(due to different effective offsets for the surface pro¬ 
file for electrostatic and Casimir interactions) may 
be avoided. This is particularly true if the plate 
is smooth enough (atomically flat) on all length 
scales so that its roughness profile does not in¬ 
fluence the distance at contact. Then, the local 
variation in the Casimir force is mainly determined 
by deviations of the surface curvature radius from 
its global value, which amounts to 1-5 %, depend¬ 
ing on the distance. These geometrical variations, 
again, are to be seen on the ‘locality scale’ of the 
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(a) 4310A (b) 4320A 

Figure 9 : (color online) Histograms of the fitted local radii on spheres of type 4310A (a) and 4320A (b), before (upper 
figures) and after (lower figures) sputtering the surfaces. In order to be comparable, the cap radii have been normalized 
by the globally fitted radius for each sphere (1 on the ordinate). The width of the histogram bars has been set to <jr/2 , 
where an is the standard deviation of the respective (assumed Gaussian) distribution. 


interaction and can only be determined from un- 
filtered AFM topography data taken at the very 
point of the force measurement. 

The present work shows that commonly used micro¬ 
spheres may feature surface corrugations which influ¬ 
ence short-scale surface interactions locally. The very 
realization of the geometry at the point of closest ap¬ 
proach determines the force quantitatively and qual¬ 
itatively. A theory which is compared to data from 
such experiments must take into account the real to¬ 
pography of the surface, as statistical measures may 
not be sufficiently accurate. 
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